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Abstract 

The search for cosmic antideuterons has been proposed as a promising method to 
indirectly detect dark matter, due to the very small background flux from spallations 
expected at the energies relevant to experiments. The antideuteron flux from dark 
matter annihilations or decays is, however, severely constrained by the non-observation 
of an excess in the antiproton-to-proton fraction measured by PAMELA. In this paper 
we calculate, for representative dark matter annihilation and decay channels, upper 
limits on the number of antideuteron events at AMS-02 and GAPS from requiring 
that the associated antiproton flux is in agreement with the PAMELA data. To this 
end, we first analyze in detail the formation of antideuterons in the coalescence model 
using an event-by-event montecarlo simulation and using data from various high energy 
experiments. We find that the resulting coalescence momentum shows a dependence 
on the underlying process and on the center of mass energy involved. Then, we 
calculate, using a diffusion model, the flux of antideuterons at the Earth from dark 
matter annihilations or decays. Our results indicate that, despite the various sources 
of uncertainty, the observation of an antideuteron flux at AMS-02 or GAPS from dark 
matter annihilations or decays will be challenging. 

Introduction 

The search for antimatter particles from dark matter annihilations or decays is usually hin- 
dered by the existence of large astrophysical backgrounds, which make difficult to inequiv- 
ocally attribute any possible excess in the fluxes to an exotic origin. A possible exception 
is the search for cosmic antideuterons. The flux of antideuterons expected from spallations 
of cosmic rays on the interstellar medium is peaked at a kinetic energy per nucleon Tg ~ 5 
GeV/n and rapidly decreases at smaller kinetic energies [TJI2]. In contrast, the spectrum of 
antideuterons from dark matter annihilations or decays is usually much flatter at low kinetic 
energies and could easily overcome the astrophysical background [3j . For this reason, it has 
been argued that the search of antideuterons from dark matter annihilations or decays is 
practically background free and that the observation of one single cosmic antideuteron would 
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constitute an evidence for an exotic antideuteron source. Several authors have calculated 
the expected antideuteron flux at the Earth in scenarios of dark matter annihilations [21 HI |5] 
or decays [6], using a factorized approximation of the coalescence model; these analyses have 
been recently revisited in [3 El E] using an event-by-event simulation. 

Antideuterons are produced by the coalescence of one antiproton and one antineutron. 
It is then apparent that there is a strong correlation between the cosmic antideuteron flux 
and the cosmic antiproton flux. In this paper we will analyze the prospects to detect 
primary antideuterons at AMS-02 [10] or GAPS [IT] , in view of the existing limits on the 
dark matter annihilation cross section or decay width into antiprotons from the PAMELA 
data on the cosmic antiproton-to-proton fraction [12]. We will show that, although the 
limits on the cross section or the decay width that follow from the non-observation on 
an excess in the measured antiproton-to-proton fraction are very sensitive to the choice of 
the propagation parameters and, to a lesser extent, to the choice of the dark matter halo 
profile, due to the correlation between the antideuteron and the antiproton production and 
propagation mechanisms, most of the uncertainties cancel out when calculating an upper 
limit on the antideuteron flux and the number of antideuteron events, thus giving fairly 
robust conclusions. 

We will analyze in section [2j the antideuteron production employing the coalescence 
model, which has only one free parameter, the coalescence momentum. We will determine 
this parameter using an event-by-event montecarlo simulation, analyzing data from various 
experiments measuring antideuteron or deuteron production. Our calculation suggests that 
the coalescence momentum is not universal, but depends on the underlying process and on 
the center of mass energy. Using the same event-by-event montecarlo, we then calculate in 
section [3] an upper limit on the antideuteron flux and the number of events at AMS-02 and 
GAPS for two representative scenarios of dark matter annihilation and decay, imposing the 
requirement that the cosmic antiproton-to-proton ratio is in agreement with the PAMELA 
measurements. Lastly, in section [4] we will present our conclusions. 



2 The coalescence model 

To describe the antideuteron production we will employ the coalescence model [UH [HJ [151 
[7J, which postulates that the probability of the formation of an antideuteron out of an 
antiproton-antineutron pair with given four-momenta kp and k£ can be approximated as 
a narrow step function 6 (A 2 + pi), where A M = kp — k£. In this model, the coalescence 
momentum p is the maximal relative momentum of the two antinucleons that still allows 
the formation of an antideuteron. 

One can show [7J that for \k D \ ^> p , being kf, = kp + k n , this ansatz leads to the 
following differential antideuteron yield in momentum space: 

&N D t 1 4 3 d 3 N p d 3 N n (k D k D \ 

75 W M = 8 ' r*» ' ^lc%a¥kr [T> T) ' (1) 

so the number of antideuterons scales with pjj. Note furthermore, that under the assumption 
of statistical independence of the antiproton and antineutron production, eq. ([lj reduces 
again to the formula used in the factorized approximation of the coalescence model [HI 
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Ill G3 E]|j In absence of a microscopic understanding of the coalescence mechanism, the 
coalescence momentum po should be determined from experiments. 

We have analyzed various high-energy processes of antideuteron production and we have 
implemented the coalescence model in an event-by-event simulation. Concretely, we have 
used the Montecarlo generators PYTHIA 6 [TB] or PYTHIA 8 [T7| to simulate the processes 
and we have selected those events with present an antiproton-antineutron pair in the final 
state with a relativistic invariant momentum difference —A 2 < p$. Furthermore, we have 
also required that the antiproton and the antineutron are close in position space, as the 
nuclear forces that are responsible for the binding of the antideuteron have a finite range 
in the order of a few fm. Hence, we have included in our analysis the antiprotons and 
antineutrons produced directly in the hadronization process as well as those produced in 
strong decays, such as from A decay which promptly decay with a lifetime ta = 5 x 10~ 24 
s very close to the primary vertex. In contrast, we have neglected those produced in weak 
decays, such as from A and £ decay, which have a lifetime ta e ~ 10 -10 s and which then 
have a macroscopic decay length, much larger than the typical range of the nuclear forces. 
In consequence, antinucleons produced in weak decays have no chance to coalesce with 
antinucleons produced in the hadronization or in other weak decays, even if their relative 
momentum is smaller than pq. 

The processes we have analyzed and the corresponding coalescence momentum inferred 
from the procedure described above are: 

• Multiplicity of antideuterons from Z boson decay measured by ALEPH [IB] . The 
number of antideuterons per hadronic Z decay with a momentum between 0.62 and 
1.03 GeV/c and a production angle satisfying \cos9\ < 0.95 was measured to be 
(5.9 ± 1.8 ± 0.5) x 10 -6 . We reproduce this number by choosing po = 192 MeV in our 
simulation (using PYTHIA 8); the quoted error can be translated to a uncertainty of 
~ 30 MeV in the coalescence momentumJ3 

• Antideuteron production in the decay of the T (IS) meson produced in e + e~ collisions, 
measured by CLEO [19]. In the data and therefore also in our simulation (using 
PYTHIA 8), the number of antideuterons is normalized to the number of T (IS) 
decaying to either ggg or ggj, denoted as 'direct' decays in [TH]. In Fig. [TJ top 
left plot, we show the data and the results of our simulation for the direct braching 
ratio into antideuterons, being the best-fit value of the coalescence momentum 
po =133 P MeV. 

• Antideuteron production in pp collisions at y/s = 53 GeV at the CERN ISR [201 12T] . 
Here the invariant differential cross section for the inclusive D production at 8cm = 90° 
was measured as a function of the transverse momentum. As shown in Fig. [T], top 
right plot, we reproduce the spectrum with PYTHIA 8 by choosing p =152 MeV. 

• Antideuteron production in e~p collisions at y/s = 318 GeV measured by the ZEUS 
collaboration [22]. The experimental analysis as well as our simulation was performed 



1 In some works the factor 1/8 appearing in eq. Q is included in the definition of the coalescence 
momentum, resulting in a value of po which is one half of the coalescence momentum defined through eq. 
Q. 

In [7J |8] a value of po — 160 MeV is extracted from the Z decay measurement. However, this value is 
obtained by allowing also the antiprotons and antineutrons originating from weak decays to coalesce. 
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Figure 1: Antideuteron/deuteron spectra measured in various experiments compared to the 
prediction of the coalescence model for a value of the coalescence momentum p = 192 MeV 
(inferred from the ALEPH measurement of the antideuteron multiplicity in Z boson decay) 
and for the best fit value for each experiment. Top left: antideuteron spectrum measured 
at CLEO from T decay, top right: antideuteron spectrum measured at CERN ISR from 
pp collisions at \/s = 53 GeV, bottom left: antideuteron spectrum measured at ZEUS from 
e~p collisions at \fs = 318 GeV, bottom right: preliminary deuteron spectrum measured by 
ALICE in pp collisions at y/s = 7 TeV. See the text for details. 



for deep inelastic scattering events with an exchanged photon virtuality Q 2 > 1 GeV 2 . 
The measured quantity is the Lorentz invariant differential antideuteron yield in the 
rapidity range \y\ < 0.4. The three data points are shown in Fig. [TJ bottom left plot, 
together with our simulation using a best-fit value of po = 236 MeV. Here we used 
PYTHIA 6, as deep inelastic scattering is not included in PYTHIA 8 yet. 

• Preliminary data on deuteron production in pp collisions at \/s = 7 TeV by ALICE 
[23] . At very high energies, the production cross sections for protons and antiprotons 
are nearly equal (see e.g. [21]); it is therefore reasonable to assume the same for the 
production cross sections of deuterons and antideuterons which allows us to extract 
information about the coalescence momentum from this measurement of deuterons. 
The double differential deuteron yield , d2 j v was measured as a function of the trans- 

J ay apT 

verse momentum in the pseudorapidity range \rj\ < 0.9 and is shown in Fig. [TJ bottom 
right plot, together with our simulation using the best-fit value of po = 230 MeV. As 
discussed in [21], PYTHIA 6 (tune D6T) reproduces the (anti)proton spectrum better 
than PYTHIA 8, therefore we used the former for our simulation. 
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Figure 2: Antideuteron injection spectra expected from dark matter annihilations 
into W + W~ (left plot) and bb (right plot) for dark matter masses moM = 100 GeV and 
tbdm = 1 TeV. The solid line corresponds to a value of the coalescence momentum po = 192 
MeV, while the shaded regions show the uncertainty band corresponding to the ranges 
162 MeV <p < 222 MeV and 133 MeV < p < 236 MeV for annihilations into W + W~ and 
bb respectively. 



We have checked, for each data set, that the Montecarlo generator used for the de- 
termination of the coalescence momentum reproduces reasonably well the corresponding 
antiproton spectrum. 

As apparent from the plots, where we also show the predicted spectra for the value of the 
coalescence momentum inferred from Z decays, po = 192 MeV, it is impossible to reproduce 
the results of all these experiments with a single coalescence momentum. Understanding 
the dependence of the coalescence momentum with the underlying hard process and with 
the center of mass energy is in our opinion an interesting problem which deserves further 
investigation, nevertheless in this paper will not pursue this direction and instead we will 
concentrate on the indirect detection of antideuterons from dark matter annihilations or 
decays, bearing in mind the above-mentioned limitation of the coalescence model to describe 
the production of antideuterons. 

We will analyze in this paper the antideuteron flux from dark matter annihilations or 
decays into W + W~ or into bb. In both cases, we will take the value of the coalescene 
momentum inferred from Z decays, po = 192 MeV. This choice is well justified for dark 
matter annihilations/decays into W + W~ , irrespectively of the dark matter mass. However, 
as argued above, the coalescence momentum for dark matter annihilations or decays directly 
into partons can differ from this value. To estimate the impact of the uncertainties in the 
coalescence model on the cosmic antideuteron fluxes, we show in Fig. |2j left plot, the range of 
antideuteron injection spectra from DM DM — > W + W~ for dark matter masses m^M — 100 
GeV and hidm = 1 TeV taking p between 162 MeV and 222 MeV, the range that results 
from the uncertainty in the determination of the coalescence momentum from the ALEPH 
data. We also show, in the right plot, the uncertainty in the antideuteron spectrum for 
DM DM — > bb, also for wdm = 100 GeV and m DM = 1 TeV, taking p between 133 MeV 
and 236 MeV, which is the range of coalescence momenta we find from our analysis of the 
experimental data on antideuteron and deuteron production. 
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3 Antideuterons from dark matter annihilations / decays 



We will assume that the Milky Way is embedded in a spherically symmetric dark matter 
halo with radial distribution p(r). In order to evaluate the uncertainity in the predictions 
of the antideuteron flux from our ignorance of the structure of the dark matter halo, we will 
show results for the Navarro- Frenk- White (NFW) density profile j25j [26] : 
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with scale radius r s = 24.42 kpc, the Einasto profile [271 1281 

Pdm(» = Poexp 

with a = 0.17 and r s = 28.44 kpc, and the much shallower isothermal profile 
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with r s = 4.38 kpc. In all the cases, the overall normalization factor po will be chosen 
to reproduce the local dark matter density p = 0.39 GeV/cm 3 [311 [32j [33J, |3U [35] with 
r Q = 8.5 kpc. 

The annihilation or decay of dark matter particles with mass wdm at the position r*with 
respect to the center of the Galaxy produces antideuterons with a rate per unit of kinetic 
energy and unit volume given by 



Qd{T,t) 



2 {(TV) f ^ (annihilations) 
" DM<r1 (decays) , 



IPdmCO 

2 m DM 



(5) 



where (cry) / denotes the velocity weighted annihilation cross-section, T f the decay rate and 
dN f D /dT the energy spectrum of antideuterons produced in the annihilation or decay channel 
/. The antideuteron injection spectra dN^/dT are calculated with PYTHIA 8 using the 
event- by-event coalescence model with po = 192 MeV, as discussed in the previous chapter. 

After being produced, antideuterons propagate through the Milky Way in a complicated 
way before reaching the Earth. We will simulate, following [SB], the antideuteron propaga- 
tion by means of a stationary two-zone diffusion model with cylindrical boundary conditions. 
Then, the number density of antideuterons per unit kinetic energy at the position r, /g(T, f), 
satisfies the following transport equation: 







df. 



D 



dt 



V • (K(T, r)Vf D ) - V ■ (V c (f)f D ) - 2hS{z) (r ann + T 



non— ann 



f B + Q D (T,f). (6) 



In this equation energy losses and reacceleration have been neglected. We will impose as 
boundary conditions that the number density of antideuterons vanishes at the boundary of 
the diffusion zone, which we approximate by a cylinder with half-height L — 1 — 15 kpc and 
radius R = 20 kpc. 
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Model 


5 


K (kpcVMyr) 


L (kpc) 


V c (km/s) 


MIN 


0.85 


0.0016 


1 


13.5 


MED 


0.70 


0.0112 


4 


12 


MAX 


0.46 


0.0765 


15 


5 



Table 1: Astrophysical parameters compatible with the B/C ratio that yield the minimal 
(MIN), median (MED) and maximal (MAX) antideuteron flux. 



The first term on the right-hand side of the transport equation is the diffusion term, 
which accounts for the propagation through the tangled Galactic magnetic field. The dif- 
fusion coefficient K(T, r) is assumed to be constant throughout the diffusion zone and is 
parametrized by: 

K(T) = K (3 H s , (7) 

where (3 = v/c and 1Z is the rigidity of the particle, which is defined as the momentum 
in GeV per unit charge, TZ = p(GeV)/Z. The normalization Kq and the spectral index 5 
of the diffusion coefficient are related to the properties of the interstellar medium and can 
be determined from the flux measurements of other cosmic ray species, mainly from the 
Boron to Carbon (B/C) ratio |37J. The ranges of the astrophysical parameters that are 
consistent with the B/C ratio and that produce the minimal (MIN), median (MED) and 
maximal (MAX) antideuteron fluxes are listed in Table [I] The second term is the convection 
term, which accounts for the drift of charged particles away from the disk, induced by the 
Milky Way's Galactic wind. It has axial direction and is also assumed to be constant inside 
the diffusion region: V c (f) = V c sign(z) k. Lastly, the third part describes annihilations 
of antideuterons with ordinary matter in the Galactic disk, which is assumed to be an 
"infinitely" thin disk with half-width h = 100 pc, and non-annihilating interactions of 
antideuterons with interstellar matter. In the latter processes antideuterons lose energy, 
leading to a redistribution of the flux towards lower energies. However, this "tertiary" 
contribution to the low energy part of the spectrum turns out to be negligibly small for 
primary antideuterons, hence we will treat these antideuterons as if they simply disappear 
from the flux. The reaction rates for these two processes can be written as 

r ann + r no „ n = (n H + 4 2 / 3 n He ) (o%? + ° D °™) v , (8) 

where un ~ 1 cm -3 is the number density of Hydrogen nuclei in the Milky Way disk and 
^He ~ 0.07 nn the number density of Helium nuclei. In this expression it has been assumed 
that the inelastic cross section between an antideuteron and a helium nucleus is related to 
the inelastic cross section between an antideuteron and a proton by the simple geometrical 
factor 4 2//3 . Following [2J, we use for our numerical analysis the estimation cr^ 1 (Tg/n) ~ 
2°|5p Cl (Tp = T D /n), where a mcl = a ann + o- non - ann and Tg/n is the kinetic energy per nucleon 
of the antideuteron. The inelastic antiproton-proton cross section is taken from [38] : 

a™ 1 (T p ) = 24.7 (1 + 0.584 Tr a115 + 0.856 Tr - 566 ) mbarn , (9) 

with Tp in units of GeV. 



7 



Finally, the interstellar flux of primary antideuterons at the Solar System from dark 
matter annihilation or decays is given by: 



where v is the antideuteron velocity. This is not, however, the antideuteron flux measured 
by experiments, which is affected by solar modulation. In the force field approximation [39] 
the effect of solar modulation can be included by applying the following simple formula 
that relates the antideuteron flux at the top of the Earth's atmosphere and the interstellar 
antideuteron flux [30] : 



where Tj$ = Ttoa + 4>f, with Tig and Ttoa being the antideuteron kinetic energies at the 
heliospheric boundary and at the top of the Earth's atmosphere, respectively, and (fip being 
the Fisk potential, which varies between 500 MV and 1.3 GV over the eleven-year solar 
cycle. For our numerical analysis we will choose 4>f = 500 MV. 

We will also make use in this paper of the limits on the dark matter annihilation cross 
sections or decay rates from the PAMELA p/p data. We have calculated these limits prop- 
agating the primary antiprotons using the same formalism as for antideuterons, with the 
obvious changes in the particle mass and a 1 ^ —> cr^f 1 (see also [4TJ EH])- 

The best present limit on the cosmic antideuteron flux was set by BESS in the range of ki- 
netic energy per nucleon 0.17 <T< 1.15GeV/n, $5 < 0.95 x 10~ 4 m~ 2 s _1 sr~ ^eV" 1 [B], 
and lies well above the expected fluxes from secondary production or from dark matter 
annihilation/decay processes. Interestingly, the sensitivity of experiments to the cosmic 
antideuteron flux will increase significantly, by more than two orders of magnitude, in the 
near future. The Alpha Magnetic Spectrometer Experiment (AMS-02) on board of the 
international space station is currently searching for cosmic antideuterons in two energy 
windows, 0.2 < T < 0.8 GeV/n and 2.2 < T < 4.4 GeV/n, with an expected flux sensi- 
tivity after five years $5 = 5 x 10 -7 m -2 s~ 1 sr~ 1 GeV~ in both energy windows [45]. Fur- 
thermore, starting in 2014, the balloon borne General Antiparticle Spectrometer (GAPS) 
will undertake a series of flights at high altitude in the Antarctica searching for cosmic 
antideuterons. In the first phase, a long duration balloon (LDB) flight will search for an- 
tideuterons in the range of kinetic energy per nucleon 0.1 < T < 0.25 GeV/n, with a 
sensitivity $5 = 6.0 x 10~ 7 m _2 s _1 sr _1 GeV~ 1 , while the second, the ultra long duration 
balloon (ULDB) flight will search in the range 0.1 < T < 0.25 GeV/n, with a sensitivity 
= 1.75 x 10~ 7 m~ 2 s -1 sr _1 GeV _1 [45J. These two experiments approach the required 
sensitivity that might allow the indirect detection of dark matter, as illustrated in Fig. [3] 
where we show the maximal antideuteron flux from dark matter annihilations into W + W~ 
and into bb which is consistent with the PAMELA measurements of the antiproton-to-proton 
fraction for m DM = 100 GeV and 1 TeV, as well as the predicted cosmic antideuteron back- 
ground from spallation of cosmic rays on the interstellar medium; for the latter, we use the 
calculation provided by [21 . The band reflects the uncertainty in the prediction due to the 
choice of the propagation model. It is important to note that, despite the predictions for the 
antiproton and the antideuteron fluxes at Earth can vary by almost two orders of magnitude 
depending on the propagation model, when normalizing the primary antideuteron flux to 
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Figure 3: Maximum antideuteron flux from dark matter annihilations into W + W~ (left 
plot) and bb (right plot) compatible with the PAMELA measurements of the antiproton- 
to-proton fraction for dark matter masses m DM = 100 GeV and m DM = 1 TeV, assuming 
a NFW dark matter halo profile and a value of the coalescence momentum p = 192 MeV. 
The red (blue) line corresponds to the MED propagation model and the red (blue) shaded 
regions spans the propagation uncertainty between MIN, MED and MAX parameters. 



the primary antiproton flux a large part of the uncertainty cancels out, due to the fact that 
the propagation of antideuterons is tightly correlated to the propagation of antiprotons. As 
a result, the uncertainty in the maximal antideuteron flux allowed by the measured p/p 
ratio amounts to less than a factor of four. It is also apparent from the plot that there is 
a strong dependence of the maximal antideuteron flux from DM DM — > W + W~ with the 
dark matter mass, in contrast to that from DM DM — > bb, due to different shape of the 
corresponing antideuteron injection spectra (c.f. Fig. [2]). 

These plots suggest that the present limits on the annihilation cross section into W + W~ 
or bb from the non-observation at PAMELA of an excess in the antiproton-to-proton fraction 
barely allows the observation of antideuterons at AMS-02 or GAPS. The interplay between 
the antiproton limits and the projected sensitivity to antideuterons at future experiments 
is further investigated in Fig. |4| where we show the dark matter annihilation cross section 
(top plots) and decay width (bottom plots) necessary for the observation of antideuterons 
originating in dark matter annihilations or decays into W + W~ (left plots) or bb (right plots) 
at AMS-02 or GAPS ULDB with a 95% confidence level (C.L.) excess over the background, 
for dark matter masses in the range 100 GeV to 2 TeV for annihilations into W + W~ 
and in the range 30 GeV to 2 TeV for annihilations into bb (the ranges for the case of 
decays are twice the ranges for annihilations). In this plot we have adopted the NFW dark 
matter halo profile, the MED propagation model and a value of the coalescence momentum 
Po = 192 MeV. In the case of annihilation or decays into W + W~ , GAPS ULDB is the most 
sensitive experiment for dark matter masses below 200 GeV, while AMS-02 is more sensitive 
for larger masses. In this annihilation channel, the antideuteron spectrum is peaked at a 
kinetic energy per nucleon T ~ mDivi/100, and rapidly decreases as the energy is reduced. 
Therefore, for large dark matter masses the number of antideuterons which fall into the 
AMS-02 high energy window 2.2 < T < 4.4 GeV/n, compensates the larger sensistiviy of 
GAPS, which is sensitive only to lower energetic antideuterons, 0.1 <T< 0.25 GeV/n. In 
contrast, the antideuteron energy spectrum for annihilation/decays into bb is much flatter, 
thus being GAPS the most sensitive experiment for the mass range of interest in this paper. 
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Figure 4: Annihilation cross section and inverse decay rate necessary for the observation 
of antideuterons at the 95% confidence level from dark matter annihilations or decays at 
AMS-02 (solid line) or at GAPS ULDB flight (dashed line) assuming the NFW halo profile, 
the MED propagation model and po = 192 MeV; the shaded region is excluded from the 
non-observation by PAMELA of an excess in the p/p fraction. The top plots correspond 
to the annihilations DM DM -> W + W~ (top left) and DM DM ->■ bb (top right), while the 
bottom plots to the decays DM — > W + W~ (bottom left) and DM — > bb (bottom right). 



We also show in the figure the 95% CL. excluded region of parameter space by the non- 
observation by PAMELA of an excess in the cosmic antiproton-to-proton fraction, adopting 
the background antiproton field derived in [36] and the proton flux measured by AMS- 
01 [37]. As apparent from the plot, the PAMELA experiment already covers, for the set of 
assumptions listed above, the whole parameter space accessible by AMS-02 and most of the 
parameter space which will be explored by GAPS ULDB. 

For the experimental searches it is more meaningful to calculate the upper limit on the 
number of events rather than the cross section. To this end, we show in Fig. [5] the maximum 
number of antideuteron events expected at AMS-02 from the 95% C.L. upper limit on the 
annihilation cross section (upper plots) or on the decay rate (lower plots) of dark matter 
particles into W + W~ and bb from the PAMELA data on the p/p ratio. We also show as 
reference the number of events expected from secondary production as well as the number 
of events necessary to claim an excess over the astrophysical background with a 95% C.L. 
Furthermore, and in order to evaluate the impact of the choice of the propagation model in 
our conclusions, we show in the plot the results for the MAX, MED and MIN propagation 
models (from top to bottom). The weak dependence of the antideuteron flux normalized to 
the antiproton-to-proton fraction becomes manifest in the plot when calculating the number 
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of events at AMS-02, and amounts to an uncertainty which is less than a 10%. It follows from 
the figure that, for dark matter annihilation/decays into W + W~, the number of expected 
antideuteron events is below one for a wide range of dark matter masses, except when the 
dark matter particle is relatively light, wdm ~ 200 GeV, where one event could be observed 
after 5 years of data taking. It is important to emphasize that the observation of one cosmic 
antideuteron at AMS-02, although a fantastic discovery in itself, cannot be interpreted 
as a signal of dark matter annihilation/decay since it could be an upper fluctuation of 
the antideuteron background. Instead, in order to claim a 95% CL. indication for an 
antideuteron excess, at least 2 events should be observed. Similar conclusions can be drawn 
when the dark matter particle annihilates or decays into bb^\ 

The prospects for detection of antideuterons from dark matter annihilations/decays at 
GAPS ULDB are shown in Fig. [6} For GAPS ULDB it is expected the observation of = 
0.04 antideuterons from spallations of cosmic rays on the interstellar medium, thus being 
practically a background free experiment. Indeed, as shown in the figure, the observation 
of one single antideuteron event would suffice to assert with a 95% C.L. the exotic origin of 
that antideuteron. Again, we find that for a wide range of parameters, it is expected less 
than one antideuteron event for the MED and MIN propagation models, except for very low 
dark matter masses. Only for the MAX propagation model, the observation of one or two 
antideuteron events would be possible in view of the constraints from PAMELA on these 
annihilation/decay channels from the non-observation of an excess in the p/p ratio, provided 
wdm < 125 GeV (m DM < 400 GeV) for annihilations into W + W~ (bb); for decaying dark 
matter the results are analogous, replacing moM with 2moM- 

It should be stressed that these conclusions apply for a value of the coalescence momen- 
tum po = 192 MeV. However, as argued in section |2j for annihilations into bb the coalescence 
momentum could be larger for heavier dark matter particles, when the center of mass energy 
of the process increases. If this is the case, the antideuteron flux, which is proportional to 
the third power of the coalescence momentum, could be significantly enhanced, improving 
the prospects to detect antideuterons from dark matter annihilations/decays at AMS-02 or 
at GAPS. For example, if the coalescence momentum for the annihilation of dark matter 
particles into bb reaches 236 MeV, as hinted by the ZEUS data, the number of expected 
antideuterons would almost double the prospects presented in the right plots of Figs. |5j [6j 

In all the plots presented in this paper we have assumed a Fisk potential <Pf = 500 MeV. 
However, data taking at future experiments will not necessarily occur during the minimum 
of the solar cycle, as was the case of the PAMELA experiment. However, increasing the 



Fisk potential in Eq. (11) has the effect of reducing the number of antideuterons at the 
top of the atmosphere, therefore the prospects presented in this paper can be regarded as 
conservative. 

Lastly, we would like to briefly discuss the dependence of the results on the choice of the 
halo model. We show in Fig|7|the maximum number of antideuteron events at GAPS ULDB 
compatible with the PAMELA antiproton-to-proton fraction for the NFW, the Einasto and 
the Isothermal profile; in each case we use the MED propagation model and po = 192 
MeV. Again, due to the fact that the primary antideuteron flux is correlated to the primary 



3 The number of background events just in the low energy window of AMS-02 amounts to 0.034, hence the 
search for dark matter particles in the window 0.2 < T < 0.8 GeV/n is background free with the AMS-02 
sensitivity. Nevertheless, restricting the search to this window yields a number of antideuterons from dark 
matter annihilations/decays which is smaller than one in the whole range of masses. 
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Figure 5: Maximal number of antideuteron events expected at AMS-02 compatible with 
the PAMELA p/p data from the annihilation channels DM DM — > W + W~ (top left) and 
DM DM — > bb (top right), and the decay channels DM — > W + W~ (bottom left) and 
DM — > bb (bottom right), assuming the NFW profile. The three dashed lines correspond 
(from top to bottom) to the MAX, MED and MIN propagation models. We also show the 
number of antideuterons expected from spallations and the total number of antideuterons 
corresponding to a 95% C.L. excess. 
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Figure 6: The same as Fig|5]but at GAPS ULDB. 



antiproton flux, the uncertainty resulting from our ignorance of the halo model almost 
exactly cancels out, giving an uncertainty in the number of antideuterons at GAPS ULDB 
which is smaller than a 50% in the case of annihilations and smaller than 5% in the case of 
decays; the uncertainty at AMS-02 is even smaller. 



4 Conclusions 

We have studied the prospects to observe antideuterons at AMS-02 and GAPS from dark 
matter annihilations or decays, in view of the stringent limits on the existence of an ex- 
otic component in the cosmic antiproton-to-proton ratio which follow from the PAMELA 
measurements. The antideuteron production was simulated using the coalescence model in 
an event-by-event montecarlo, being the coalescence momentum the only free parameter. 
To determine the value of this parameter we analyzed data from the experiments ALEPH, 
CLEO, CERN-ISR and ZEUS, as well as preliminary data from ALICE. We have found 
that the event-by-event montecarlo reproduces fairly well the data, however the coalescence 
momentum is not a universal parameter but depends on the underlying process and the 
center of mass energy. This result suggests that the coalescence model commonly used in 
the literature is oversimplified and should be extended. In this paper we have not attempted 
to construct an improved coalescence model but, in view of the qualitatively good descrip- 
tion of the experimental data on antideuterons, we have used the simplest model adopting 
the value of the coalescence momentum inferred from the ALEPH measurements of the 
antideuteron multiplicity in Z boson decay and we have discussed the impact of a different 
value of the coalescence momentum on our conclusions. 
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Figure 7: Sensitivity to the choice of halo profile of the maximal number of antideuterons 
expected from dark matter annihilations (top plots) or decays (bottom plots) into W + W~ 
(left plots) or bb (right plots) at GAPS ULDB. The three lines correspond (from top to 
down) to the Isothermal, the NFW and the Einasto profile. 
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The calculation of the maximal antideuteron flux and the maximal number of an- 
tideuteron events from dark matter annihilation or decays at AMS-02 or GAPS suffers 
from various sources of uncertainty. To evaluate the particle physics uncertainties we have 
considered two representative dark matter annhihilation and decay modes, namely into 
W + W~ and into bb. Furthermore, we have taken into account various astrophysical uncer- 
tainties, such as the choice of the dark matter halo model, of the propagation model or the 
value of the Fisk potential. While the primary antiproton and antideuteron fluxes are very 
sensitive to the choice of the propagation model and, to a lesser extent, to the choice of the 
halo model, due to the fact that the antiproton and antideuteron propagation are tightly 
correlated, most of the uncertainties cancel out when calculating the primary antideuteron 
flux normalized to the primary contribution to the cosmic antiproton-to-proton ratio, giving 
maximum fluxes which are fairly insensitive to the astrophysical uncertainties. 

We have found that the stringent limits on the abundance of primary antiprotons in 
cosmic rays from the PAMELA experiment severely constrain the possibility of observing 
antideuterons from dark matter annihilations or decays at AMS-02 or GAPS. At AMS-02 the 
expected number of antideuterons from dark matter annihilation or decays is smaller than 
one, except when the dark matter particle is relatively light, where at most one event could be 
found, according to our calculations. Furthermore, at AMS-02 we expect the observation of 
~ 0.2 secondary antideuterons, hence even if the upper limit on the number of antideuterons 
is saturated, the observation of one single antideuteron could not be inequivocally attributed 
to dark matter annihilations or decays. The conclusions for the GAPS experiment are 
analogous for the MIN and MED propagation model and more optimistic for the MAX 
propagation model. 
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